Odd parity charge density-wave scattering in cuprates 
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We investigate a model where superconducting electrons are coupled to a frequency dependent 
charge-density wave (CDW) order parameter A r (cj). Our approach can reconcile the simultaneous 
existence of low energy Bogoljubov quasiparticles and high energy electronic order as observed in 
scanning tunneling microscopy (STM) experiments. The theory accounts for the contrast reversal 
in the STM spectra between positive and negative bias observed above the pairing gap. An intrinsic 
relation between scattering rate and inhomogeneities follows naturally. 



PACS numbers: 71.10.Hf, 74.72.-h, 74.25. Jb 



The phenomenon of charge ordering and its relation to 
superconductivity remains a puzzling issue in the physics 
of iiigh-T c cuprates. Whereas the occurence of charge in- 
homogeneities was early on evidenced by local probes and 
discussed in the context of electronic phase separation [l| 
Tranquada and coworkers Q succeeded in probing a one- 
dimensional static spin and lattice modulation in a rare 
earth codoped lanthanum cuprate compound (LCO) by 
elastic neutron scattering. The unambigous existence of 
an associated charge modulation has been established 
only more recently with resonant soft X-ray scattering^]. 
In other cuprate compounds the existence of charge or- 
dering so far comes from surface sensitive probes, like 
scanning tunneling microscopy (STM) and angle-resolved 
photoemission spectroscopy (ARPES). STM experiments 
performed on bismuthate and oxychloride superconduc- 
tors see a complex modulation of the local density of 
states (LDOS) both in the superconducting (SC) state 

0, 0, |E 0, H, M and above T c B ES ELU- In both cases 
one observes peaks in the Fourier transform of the real 
space LDOS at wave-vectors Q — 2n / (Aao) . . .2tt / (5ao) 
suggestive of checkerboard or stripe charge order. How- 
ever, the debate is about the question whether these 
peaks are non-dispersive in energy (and thus signature 
of 'real' charge order) or follow a bias-dependent disper- 
sion due to quasiparticle interference (QPI). In the latter 
case the spatial LDOS variations can be understood from 
the so-called octet model 0, EH which attributes the 
modulations to the elastic scattering between the high 
density regions of the Bogoljubov 'bananas' in the su- 
perconducting state. Recent STM investigations IS 2A\ 
may resolve this apparent conflict since they suggest that 
both, dispersive and non-dispersive scattering originates 
from different regions in momentum and energy space. 
The states in the nodal region which are well defined in 
k-space and undergo a transition to a d-wave SC state 
below T c are then responsible for the low energy QPI 
structure of the LDOS, whereas the ill-defined k-space 
'quasiparticle' states in the antinodal regions are respon- 
sible for the non-dispersive charge order above some en- 



ergy scale fio- It is therefore a key issue to understand 
the nature of this charge ordering, which seems particu- 
larly elusive, at least at low energies, both in STM and 
ARPES experiments. 

In this paper we propose a phenomenological model 
which captures the above physical scenario by consid- 
ering a frequency dependent charge-density wave order 
parameter analogous to the frequency dependent super- 
conducting order parameter of Eliashberg theory. By 
construction the order parameter vanishes on the Fermi 
surface so that for small energies the system appears ho- 
mogeneous and the concept of QPI applies while for large 
energies it appears to be charge ordered. As a bonus 
the system shows a strong high energy contrast reversal 
of the LDOS as a function of energy in agreement with 
experiments 24j and an intimate relation between inho- 
mogeneities and quasiparticle scattering rate. 

We consider here a two-dimensional system of SC itin- 
erant electrons scattering with a charge density wave 
order parameter with an internal dynamics. This is 
phenomenologically represented by a retarded local self- 
energy where the space and frequency dependencies are 
factorized 



-.CDW 



(co)=A r (uj) = A Q r +v 2 J(Lu). 



(1) 



A conventional CDW without internal dynamics has 
f(ui) — and A° ^ 0. The function f(u>) describes 
the internal dynamics of the CDW, that is particle-hole 
correlations building the charge density wave are time 
dependent. The spatial dependence is implemented via 
the modulation of v 2 . We will restrict to Im/(u;) > 0, 
and v r real, which ensures ImS^ DW > as required for 
stability. The internal dynamics can be motivated on the 
basis of a frequency dependent interaction, presumably of 
electronic origin and the same that originates supercon- 
ductivity but which we do not need to specify. It implies 
that charge piling is retarded which is physically appeal- 
ing since it reduces the Coulomb penalty as in Eliashberg 
theory of superconductivity. Such retardation effects are 
especially important for cuprates which are characterized 



2 



by a large Coulomb repulsion. 

In order to implement correctly the analytical proper- 
ties of the self-energy (Kramers-Kronig, etc.) it is conve- 
nient to make a pole expansion, 



/H = £- 



i 



This maps the self-energy to that of an effective Fano- 
Anderson model of itinerant electrons with conduction 
bandwidth of order t which can locally hop onto a distri- 
bution of / levels via a sitc-dcpcndcnt hybridization term 
v r . The bath of / levels simulates the scattering due to 
the charge fluctuations and provides a Hamiltonian for- 
mulation for the self-energy. 

We will consider the dynamic case with A" = and 
compare with the conventional CDW model (v T = 0, 
A" ^ 0). A small A^! component added to the fre- 
quency dependent case will not change our results sig- 
nificantly. In addition we will take the modulation to 
be periodic vl — 1/N S J2n( v2 )nQ exp(inQr). Just as a 
superconductor with a frequency dependent order pa- 
rameter has long range off-diagonal order, even with 
A£ = 0, the dynamical model has diagonal order given by 
n r = l/N s Y^, nt k J dwImGk,k+nQ cos(nQr). In real ma- 
terials the charge modulations will be linearly coupled to 
the impurity potentials and the system will easily lose 
long range order but keeping short range modulations 
with translational symmetry breaking. This important 
effect can be incorporated but will be neglected for sim- 
plicity. 



Motivated by the experiments by Kohsaka et al. [15 1 
and previous evidence on stripes [2| we consider a one- 
dimensional modulation with wave- vector |Q| = 27r/4 
along the x— direction. Results for checkerboard patterns 
will be similar except for the absence of C4 symmetry 
breaking. The upper right insets to Fig. [3^,,b display the 
associated modulation A° and v%- The resulting charge 
modulation can be seen as bond-centered (hole) stripes 
separated by 4 lattice constants respectively. Our main 
conclusions do not depend on this particular choice. 

The Green function (GF) in k-space can be obtained 
from the coupled system of equations 



[w - £k+nQ - «o/M] Gk+nQ,k+mQ - S nm 

,k+mQ (2) 

where the dispersion £k is measured with respect to the 
chemical potential. While the zero momentum compo- 
nent of the static CDW order parameter A° =0 gets triv- 
ially reabsorbed in the chemical potential the same does 
not occur with the dynamical part. That is T,q DW (uj) = 
VqJ(uj), with Vq = (w 2 ) q= o, has to be explicitely taken 
into account as can be seen from the mapping to the 
hamiltonian structure which guarantees the preservation 



of sum rules for the spectral function. Furthermore, as 
discussed above, we are forced to take v% > which im- 
plies that the Fourier component Vq has to be positive for 
any non-zero modulation. Thus we have the surprising 
result that the dynamical CDW scattering produces a 
momentum independent scattering rate. For a given am- 
plitude of modulation, there is a lower bound for such a 
scattering rate which is determined by taking the smallest 
v% to be zero. This leads to a scattering rate of the charge 
carriers by the CDW fluctuations T° DW ~ Im^ flly (w). 
In the following we make this minimum choice which 
makes r^ DW to be explicitely determined by the sum of 



amplitudes of the electronic inhomogeneity 



Such 



an intrinsic relation between electronic inhomogeneity 
and inelastic scattering rate has been recently revealed by 
STM experiments on Bi2212 materials [3] where it has 
been shown that the LDOS spectra can be parametrized 
based on a model with SC d-wave order supplemented 
by an energy dependent scattering rate r^ DOS = auu. 
The parameter a varies spatially and in the regions with 
pronounced charge order acquires values up to a w 0.4. 
Also ARPES spectra [H, [lj} from Bi2212 materials sug- 
gest a marginal Fermi liquid (MFL) type self-energy (i.e. 



ImS ~ t ar p es 



which is increasing with underdop- 



ing similar to T^ DOS . Motivated by these experimental 
findings we describe the CDW dynamics (and therefore 
T C , DW ) by a marginal Fermi liquid (MFL) type structure 

m 



f(ui) = 2wln- 



iirT 



(3) 



where T — kT corresponds to a temperature scale and 
fl denotes an upper cutoff of the boson spectrum from 
which the MFL self-energy is derived The ampli- 
tudes of the charge modulation are chosen such that 
vl = QMeV which yields ImS(w) = uglm/(w) w 0.25|w| 
at small frequencies which is close to the average r^ DOS 
observed. The dispersion £k includes hopping to near- 
est (~ t\ = 250meV) and next nearest (t% — — OAti) 
neighbors as appropriate for Bi2212 materials. 

Fig. [1^, demonstrates that for the dynamic case quasi- 
particles around the Fermi energy are protected from 
the CDW scattering. In contrast, static CDW order (cf. 
Fig. [Tb) produces the usual FS reconstruction, i.e. flat- 
tening of the dispersion around (n, 0) and gap formation 
between the nodes and (0, 7r) due to nesting. At higher 
energies the effect of the dynamic scattering on the spec- 
tral function is similar to that of a static CDW, with an 
additional broadening due to the imaginary part of the 
self-energy which is increasing with frequency. 

As mentioned above the interpretation of STM exper- 
iments, which are usually taken at very low tempera- 
tures, requires the implementation of superconductivity 
into the formalism which in Eq. [2]can be added as a BCS- 
type SC self-energy ££ c = (A| c ) 2 / [to + e k + f* (-u)Ho 
the bracket on the l.h.s. Recent STM experiments [20( 
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FIG. 1: Density plots of the integrated spectral weight in a 
window of ±5meV around Ef- (a) dynamic CDW (A? = 0); 
(b) static CDW with A° - A? = 0.054eV (cf. upper right 
inset to Fig. [3j|. Residual parameters as given in the text. 
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FIG. 2: (Color online) Full line: Anharmonic SC gap struc- 
ture Eq. Q implemented in the calculation of the LDOS; 
Dashed: Pure d-wave contribution to the anharmonic gap 
around the nodes (<&/ir = 0.25). 



have revealed that the shape of the tunneling spectra of 
Bi2212 can be fit by including higher harmonic contribu- 
tions in addition to the simple d-wave form. Following 
this observation we take 

A s h c = A-o(n) [cos(nfc x ) - cosfnfej,)] /2 (4) 



and restrict to the first three harmonics with Ao(l) = 
AOmeV, A (2) = -WmeV, and A (3) = 5meV. For 
simplicity we take this gap function to be frequency in- 
dependent. The resulting gap structure is shown in Fig. [2] 
in comparison with a pure d-wave form (Ao(l) = 18meV) 
which fits the anharmonic gap in the nodal region. The 
larger gap in the antinodes is supposed to be due to local 
pairing and persist above T c producing the pseudogap. 

In Fig. [3] we report the LDOS structure for the model 
with static (a) and dynamic (b) order parameter includ- 
ing the SC gap. The static order parameter was chosen 
in order to give a similar magnitude of the intensity dif- 
ference among different sites (contrast) at high positive 
bias as the dynamical model. 

Due to particle hole mixing superconductivity tends 
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FIG. 3: (Color online) LDOS for the model with static CDW 
scattering A§ — A? = 0.054eV (a) and frequency dependent 
CDW scattering f(u) (b). Upper insets: Sketch of the under- 
lying modulation of v^. Lower left insets: Fourier tranformed 
of the LDOS at the CDW scattering vector Q = 2£. The 
phase has been chosen such that Im/9(Q) = 0. Lower right 
insets: LDOS at u = +5meV (squares) and ui = — 5meV (cir- 
cles). Further parameters: Chemical potential fj, = —0.23eV 
(doping x « 0.07), T = lmeV and fi = leV. 



to soft the charge order related LDOS modulations at 
low energies, however the inset in (a) reveals that in the 
static case substantial residual modulations persist at an 
energy of w = ±5meV specially for negative energies. 
This contrast is strongly suppressed for the dynamic or- 
der parameter [see inset in (b)] as seen experimentally. 
This is easy to understand due to the structure of f(u>) 
which vanishes at low energies. However, this effect is 
not very apparent and would hardly allow one to decide 
for the dynamical model. A much more dramatic differ- 
ence is instead provided by the sign of the contrast when 
the bias is reversed. For the static order parameter we 
observe contrast reversal below the energy of the super- 
conducting gap (i.e. sites with high intensity and low 
intensity LDOS are inverted when one changes the sign 
of the bias) while the ordering of the intensities tends 
to be substantially preserved at high energies (contrast 
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preservation). For the dynamical case the situation is the 
opposite: there one finds contrast preservation inside and 
contrast reversal above the superconducting gap. Thus 
the dynamical model agrees with the ubiquitous contrast 
reversal found in STM at large energy concomitant with 
the absence at low energies [24] . 

How general is this result? In the static case contrast 
reversal at high energy can be achieved for special sit- 
uations (e.g. half-filled system with square Fermi sur- 
face, i.e. complete nesting) where asymmetry occurs for 
higher energies. In these cases the imaginary part of the 
off-diagonal GF IulG^+q has a sign change at the chem- 
ical potential since fc-states with ek — fJ- < only connect 
to k + Q-states above the chemical potential and vice 
versa. In case of the dynamical order parameter the sign 
change has its origin in the dynamical self-energy since 
Gk,k+Q ~ /(w) and Ref(uj) = —Ref(—u) and thus, it 
does not require special conditions as nesting. Indeed 
we have found similar results for other periodicities and 
fillings. 

The hump feature on the positive frequency side of the 
LDOS in Fig. [3] originates from the nesting of Q between 
the antinodal FS segments. This induces a concomitant 
increase of the LDOS modulation. which for our under- 
doped system occurs at energies ~ 90 meV above Ep. 

To analyze the difference between frequency dependent 
and frequency independent CDW scattering in more de- 
tail we show in the lower left insets to Figs. [3]the Fourier 
transformed LDOS at multiples at the CDW wave vector. 

As pointed out in Refs. [2ll . |22| the symmetry of the 
wave-function in a SC with static CDW order implies 
p(Q, uj) — — p(Q, — oS) for sufficiently small lu. In fact, we 
clearly observe this asymmetry in the LDOS in Fig. [3^a) 
(black, solid curve in the lower left inset) on the en- 
ergy scale of the SC gap in contrast to the experimen- 
tal observation [B| of a symmetric p(Q,w) in this regime. 
This was explained as due to a modulation of the SC 
pair density H, 23| ■ Fig. [3](b) shows that the symmetry 
is also achieved for a dynamic CDW order parameter. 
Notice that this symmetry can not persist up to high en- 
ergies since contrast reversal implies that at high energy 
the dominant Fourier component should change sign. In- 
deed our dynamic model with A CDW (u;) ~ f{u>) with 
Ref(ui) = —Ref(-oj) reconciles both aspects, low en- 
ergy symmetry of p(Q,uj) and a strong high energy odd 
component leading to contrast reversal. 

The present model applies to cuprates where the sys- 
tem at the microscopic scale shows breaking of trans- 
lational symmetry but sustains low energy quasiparti- 
cles weakly affected by that. A prominent example 
is Lai.875Bao.i25Cu04 where a surprisingly well defined 
Fermi surface has recently been observed [25] . In addi- 
tion, there are situations with no apparent translational 
symmetry breaking but with a high energy response re- 
sembling that of an ordered system. We have recently 
presented a scenario for the spectral function in this dis- 



tinct physical situation 26]. 

Summarising we have presented a simple phenomeno- 
logical model of a frequency dependent CDW order pa- 
rameter in cuprates which accounts for the dichotomy 
between low energy Bogoljubov quasiparticles and high 
energy electronically ordered states [13J with the con- 
comitant crossover from LDOS modulations determined 
from QPI to dynamical CDW scattering at large frequen- 
cies. For simplicity we have taken into account a charge 
order parameter but we expect similar results will apply 
for a magnetic order parameter. The model has a series of 
implications which agree with experimental observation: 
i) Intimate relation between quasiparticle scattering rate 
and amplitude of charge order, ii) Suppression of con- 
trast in the LDOS at low energies. Hi) Contrast reversal 
(preservation) at high (low) energies, iv) Even symmetry 
of the low energy Fourier transformed LDOS. v) Feature- 
less Fermi surface. While a naive CDW order parameter 
will be in contradiction with most of these observations 
the dynamical model reconciles the long suspected charge 
ordering with these apparently contrasting experiments. 
Our results for the scattering time suggest to search for 
the origin of marginal Fermi liquid behavior in an Eliash- 
berg treatment of charge (and) or spin density waves [27| . 
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